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Abstract 

We present the results of magnetization and magneto-transport measurements in the super¬ 
conducting state of an as-cast Nb 75 Zr 25 alloy. We also report the careful investigation of the 
microstructure of our sample at various length scales by using optical, scanning electron and trans¬ 
mission electron microscopies. The information of microstructure is used to understand the flux 
pinning properties in the superconducting state within the framework of collective pinning. The 
magneto-transport measurements show a non-Arrhenius behaviour of the temperature and field 
dependent resistivity in the flux flow region. This non-Arrhenius behaviour is understood in terms 
of a model, which was originally proposed for viscous flow of disordered solids and is popularly 
known in the literature as the ‘shoving’ model. The activation energy for flux flow is obtained from 
magneto-transport measurements and is assumed to be mainly the elastic energy stored in the 
flux-line lattice. The critical current density estimated from magnetization measurements is mod¬ 
erately high of the order of 10® Am“^ at a temperature of 2K. The scaling of pinning force density 
with respect to reduced field indicates the presence of two pinning mechanisms of different origins. 
The elastic constants of the flux-line lattice are estimated from magnetization measurements and 
are used to estimate the length scale of vortex lattice movement, or the volume displaced by the 
flux-line lattice, during flux flow. It appears that the vortex lattice displacement estimated from 
elastic energy considerations is of the same order of magnitude as that of the flux-bundle hopping 
length when a finite resistance appears during flux flow. Our results could provide possible direc¬ 
tions for establishing a framework where vortex matter and glass forming liquids or amorphous 
solids can be treated in a similar manner for understanding the phenomenon of viscous flow in 
disordered solids or more generally the pinning and depinning properties of elastic manifolds in 
random media. 

PACS numbers: 74.25.Uv, 74.25.Wx 
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I. INTRODUCTION 


The phenomena of flux creepi and thermally activated flux flow^ were initially discovered 
in conventional hard type-II superconductors. The “irreversibility line” in the H — T plane, 
which separates the reversible (or flux flow) and irreversible magnetization (or flux pinned) 
regions, was thought to be exclusive to the high-temperature superconductors and was thus 
interpreted as a “quasi de Almeida-Thouless line” in analogy to the irreversibility observed 
between the zero-field-cooled (ZFC) and field-cooled (FC) magnetization in spin glasses 
and the possibility of a superconducting glass phase was conjectured.- The vortex-glass 
phase was suggested to be a new thermodynamic phase within the mixed state of a type-II 
superconductor. - The term vortex glass was used to the signify the destruction of long range 
translational order in the Abrikosov lattice due to flux pinning by the underlying disorder 
in the superconducting material. The irreversibility line was shown to be a sharp boundary 
across a thermodynamic phase transition between the disordered vortex solid (or glass) and 
the vortex liquid where thermally activated flux flow is possible.- An alternate interpretation 
of the irreversibility line, without invoking the glass phase in the mixed state, was later 
offered in terms of the more conventional phenomenon like the flux creep. - These paradigms 
of high-Tc superconductors were applied to conventional type-II materials like niobium thin 
hlms where it was found that the irreversibility line represents the vortex-melting transition.- 
The model of the “giant flux creep”- was extended to explain the width and shape of the 
resistive transition as a function of applied magnetic field.- The key ingredient to the flux 
creep picture is the estimate of the activation energy Uq which must be surmounted for the 
flux movement to occur.- The elemental volume of the flux line lattice which moves during 
the creep is governed by the activation energy and thus gives the temperature width of the 
resistive transition as a function of applied field.-- The superconducting glass picture^ and the 
giant flux creep picture^ were thus shown to be qualitatively similar by accounting for the held 
dependence of width of transition (and in-turn the shape of irreversibility line) without the 
need for invoking sample inhomogeneities.- The combination of quenched disorder, higher 
temperatures, shorter coherence length and large magnetic penetration depth in case of 
high Tc superconductors lead to a modihcation of the mean held phase diagram of the 
conventional low Tc superconductors.- 

An alternate way of studying the vortex matter is to consider the collection of vortices 


3 


as an elastic object,- instead of treating the problem within the framework of microscopic 
Ginzbnrg-Landan phenomenology.- The rich variety of vortex phases within the snpercon- 
dncting state could be well explained by an elastic theory of vortex matter.— Ti By sup¬ 
plementing the traditional theory of collective pinningf^ with the concepts of elastic media 
in random potential, the various properties of the vortex matter could be described in an 
unihed manner. The application of the theory of elasticity to vortex matter has found simi¬ 
larities in case of other systems like magnetic domain walls, charge density waves and wigner 
crystals.- The vortex matter has thus provided a model system to study the properties of 
moving glasses or driven elastic media in presence of quenched disorder.— 

In the work presented here, we experimentally attempt to hnd similarities between the 
flux flow in vortex matter with one more class of a driven elastic system, namely the flow 
of a molecular glass or viscous liquid. We try to address one of the fundamental questions 
arising in the studies of vortex matter, which is, “How much distance does the fluxon (or the 
elemental volume of flux-line lattice) move for a hnite resistance to appear across the resistive 
transition in a hard type-II superconductor?” While this question can be of technological 
importance for tuning the nature of quenched disorder to achieve desirable critical current 
densities, it can also lead to fundamental understanding of elastic systems influenced by 
disorder and driven by external force. We have chosen the Nb75Zr25 alloy for this study. 
The Nb-Zr alloy system had generated considerable interest in the past for its perceived use 
in superconducting wires for high held magnets.— The (3 phase alloy with the b.c.c. structure 
was well studied due to its high critical current carrying ability.— The same composition 
in this alloy system was used to investigate fundamental phenomenon like the hux creep in 
hard superconductors^ or the scaling laws for hux pinning in hard type-II superconductors,— 
which are now known popularly as Kramer scaling. 

Our approach to solve the above posed question is as follows. We hrst characterize 
the microsctructure at various length scales of the Nb75Zr25 superconducting alloy under 
consideration to know about the nature of disorder which can act as pinning potential for 
the hux-line lattice. The estimate of the activation energy required for hux-line movement 
is obtained from electrical resistivity measurements across the transition as a function of 
both temperature and magnetic held. The non-Arrhenius shape of the resistive transition 
draws our attention to a model of viscous how which has been used earlier to explain a 
non-Arrhenius temperature dependent viscosity across a molecular glass transition. The 
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model is based on purely elastic energy considerations where the activation energy is related 
to the shear modulus through some sort of a correlation volume and is popularly known as 
the shoving model.— The held-dependent pinning force in the mixed state of our sample can 
be explained only if more than one pinning mechanisms are considered. Such a possibility 
is explored within the framework of collective pinning,— where the estimates of correlation 
lengths over which the hux bundles are pinned can be obtained. The elastic constants 
of vortex matter like the shear modulus (Cge) and the tilt modulus (C44) are obtained 
from magnetization measurements using this framework. The activation energy for hux- 
line movement and the elastic constants of the vortex matter are then correlated through 
the correlation volume of hux bundle to get the estimate of the distance over which the 
hux-line lattice moves before a hnite resistance appears across the transition. Interestingly, 
we hnd that the displacement of vortex matter estimated from the shoving model is of the 
same order of magnitude as that of the hux-bundle hopping length when a hnite resistance 
appears across the resistive transition. The viscous how of disordered solids and the hux-line 
movement in vortex matter appear to be just diherent cases of the same general phenomenon. 
These results could possibly lead to newer directions for encompassing diverse phenomenon 
arising in periodic systems inhuenced by quenched disorder with vortex matter as a model 
system. 

II. EXPERIMENTAL 

NbrsZrss alloy ( Igm mass) was prepared by arc melting the constituent elements of 
99.99% purity in a water cooled copper hearth placed in an inert argon atmosphere. The re¬ 
sulting sample button was remelted six times to ensure homogeneity and was not subjected to 
any further heat treatment. The sample was then characterized by x-ray dihraction (XRD) 
using a commercial dihractometer (Panalytical X-Pert PRO MRD) with Cu-Kq, radiation. 
9-29 scans were recorded in the Bragg-Brentano geometry. A monochromatized (Cu-Kqi) 
and collimated (about 20 arc-sec in the plane of scattering) x-ray beam was obtained using 
a hybrid monochromator. The obtained x-ray peak widths, which were signihcantly larger 
than the instrumental broadening, were used for the estimation of particle size and lattice 
strain. A small piece of sample was used for optical metallography which was subjected to 
slow grinding using SiC paper and polishing using diamond compound on a soft cloth. After 
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obtaining a mirror qnality finish on the snrface, the sample was chemically etched using a 
mixture of 70 ml CsHgOa (lactic acid) with 30 ml HNO 3 (nitric acid) and 2 ml HF (hy¬ 
drofluoric acid) for about 15 seconds. The microstructure was observed with a commercial 
inverted metallurgical optical microscope (Leica DMI5000M). A commercial scanning elec¬ 
tron microscope (Philips, XL30CP) equipped with a energy dispersive spectrometer (Bruker, 
XFlash@ Silicon Drift Detector) was used for observing the microstructure on sub-micron 
length scales and determining the chemical composition at various locations of the sam¬ 
ple. For visualizing the disorder/grain structure at even smaller length scales, a small piece 
of the sample was cut from the same parent button for transmission electron microscopy 
(TEM) studies. A commercial TEM (Phillips, CM200) with a tungsten filament as cathode, 
was used at an accelerating anode voltage of 200 kV. Electrical resistivity was measured 
using the standard four-probe technique with a home-made variable temperature insert in a 
commercial superconducting magnet and cryostat system (American Magnetics Inc., USA) 
Magnetization (M) measurements were performed as a function of field {H) and tempera¬ 
ture (T) using a commercial vibrating sample magnetometer (VSM Quantum Design, USA). 
We consistently use the SI units for all the physical properties reported here. 


III. RESULTS AND DISCUSSION 


Microstructure and nature of disorder 


Figure 1(a) shows the XRD pattern of the sample. The peaks which are quite broad in 
nature could be indexed with the /9 phase b.c.c. structure with a lattice constant of 3.368 
A. To estimate the grain size and lattice strain, we use the Williamson-Hall (W-H) plot,— 


which is shown in the inset of fig. 1(a) The plot is between AiF = (cos( 6 ')A( 6 '))/A and 


K = sin{6)/X, where A{6) is the full width at half maximum for each of the Bragg peaks. 
The (110) peak appears to be anomalously sharp and does not belong to the straight line 
W-H fit obtained from the other peaks. The reciprocal of the y-intersept gives the average 
grain size which is estimated to be around 98 nm. The slope of the W-H plot gives the 


lattice strain which is found to be about 1.33%. The microstructure (see fig. 1(b)) shows a 
dendritic growth during the solidification of the alloy from melt, which is quite commonly 
seen in many materials when a thermal gradient is present between the melt and the solidified 
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(a) (b) 

FIG. 1: (a) The x-ray diffraction pattern of Nb 75 Zr 25 as-cast alloy. Inset shows the Williamson- 
Hall plot, (b) The microstructure of the same sample observed under an optical microscope. The 
dendritic pattern is clearly visible. 


portion of the material.— The combined resnlts of XRD and optical metallography show that 
the alloy has probably solidihed with snbstantial disorder dne to the rapid cooling it had 
to experience because of the small size of the sample button placed in the water cooled 
copper hearth. The rapid cooling probably did not allow sufficient time for the phase to 
grow homogeneously. 

The composition variation across this dendritic pattern was determined using energy 
dispersive spectroscopy. Figure |2] shows the microstructure at two different length scales at 
a random location on the sample surface. When the composition is determined over a larger 


area (£g. 2(a)), it is quite close to the nominal composition of Nb 75 Zr 25 within the error bar 
of EDS measurements, which is nearly 3% for Nb and 1% for Zr. At certain locations on the 


protrusions of the dendritic arms (£g. 2(b)) the Zr concentration is much higher than the 
target composition (25%) and has reached up to 40 % at certain locations of the sample. 
These results imply that the sample is a mixture of Zr-rich and Zr-de£cient regions on the 
local scale even though the 3:1 ratio of Nb to Zr is maintained on a larger length scale. 
Figure [3] shows the results of the high resolution TEM studies which bring out the nature 


of disorder in a more clear manner. In fig. 3(a) the major defects are lattice plane bending, 
mismatch in interplanar spacing and edge dislocations. The inset shows the selected area 
diffraction (SAD) pattern which shows the elongation of Lane spots. It is known that the 
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(a) (b) 


FIG. 2: Scanning electron microscopy images at two different magnifications, (a) The magnifica¬ 
tion is 200x and the average composition of Nb 75 Zr 25 is observed over the large area marked as a 
rectangle, (b) The magnification is 2000x and a region on a dendritic arm is in focus. The area is 
Zr-rich in composition. 


lattice constant in the Nb-Zr alloy system varies with Zr concentration.^^ Thus the mismatch 
in interplanar spacing from 2.61 A to 2.42 A at a particular location shown in fig. 3(a 


IS 


probably due to composition variation across that region. This observation is consistent 
with the results of EDS measurements mentioned earlier. The formation of dislocations 
(shown in encircled region) is most probably due to a large lattice mismatch at certain 
locations due to compositional variation. When the strain increases beyond a certain limit 
in heterostructures, the interface accommodates a part of the lattice mismatch through 


the introduction of dislocations.-^ Figure 3(b) shows another location of the sample which 
contains lattice plane bending and possibly a jog at one of the edge dislocation. Figure 


3(c) shows a location of the sample within the dendritic growth where the lattice is heavily 


deformed due to the composition variation. It should be noted that no well defined grains 
(or grain boundaries) were observed on length scales estimated from x-ray diffraction. The 
sample is mostly a disordered alloy with intermittent disturbances in the periodicity of lattice 
at very short length scales without any sharp grain boundaries. As we see later, the typical 
length scale of these disturbances in the periodicity of lattice is smaller than the coherence 
length, which makes the collective pinning of the flux line lattice^ more probable than single 
particle pinning. The applicability of the collective pinning theory will also be shown during 
the course of this article by estimating certain superconducting parameters. 













FIG. 3: High resolution transmission electron microscopy images at different locations of the 
sample, (a) Location where the lattice mismatch and dislocations can be seen, (b) Arrows mark 
the plane bending and possible jog (middle arrow). Circles mark the edge dislocations, (c) The 
protrusions of the dendritic growth are amorphous regions that lead to plane bending. 

B. Non-Arrhenius resistive transition and activation energy 


The microstructure of the sample at various length scales shows that the material is 
highly disordered and this inhomogeneity is expected to manifest in the superconducting 
properties as well. It is known that the Tc in Nb-Zr system varies quite substantially 
with Zr concentration.^! Moreover the lack of long range ordering in the underlying lattice 
could also result in a response of the sample which is quite different from the bulk of same 
alloy system. Contrary to this expectation in the case of the present sample, we see that 
the material behaves like a typical bulk type-II hard superconductor. We highlight this 


observation before discussing the results of resistivity measurements. Figure 4(a) shows 
the temperature dependent magnetization of the Nb 75 Zr 25 sample in a low applied magnetic 
held of = 2mT. The measurement was performed under three different thermomagnetic 
histories. 

In the zero-held-cooled (ZFC) protocol, the sample is hrst cooled in zero held down 
to the lowest temperature. The magnetic held is applied at the lowest temperature and 
magnetization is measured during warming of the sample. The superconducting transition 
temperature T^ is estimated to be slightly below IIK. In the held-cooled-cooling (FCC) 
protocol, the sample is cooled in presence of held from above Tc. The temperature at which 
the FCC curve deviates from the ZFC curve is termed as the irreversibility temperature Tirr- 
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(a) (b) 

FIG. 4: (a) Low field magnetization as a function of temperature in three different thermomagnetic 
histories, (b) Electrical resistivity as a function of normalized temperature in presence of magnetic 
field. The lines are a guide to the eye. The width of the transition as a function of magnetic field 
is shown as the inset. The solid curve in the inset is a fit of experimental data with equation [TJ 


After reaching the lowest temperature in presence of held, the measurements are performed 
while warming the sample in held and this protocol is named as the held-cooled-warming 
(FCW) protocol. The FCW curve deviates from the FCC curve and meets the ZFC curve at 
a temperature which is slightly lower than These thermal history ehects, especially 
the hysteresis between the FCC and FCW curves, are indications that the sample behaves 
like a typical bulk hard type-II superconductor.— >2^ We had shown earlier that in case of 
nanocrystalline superconductors,— the FCC and FCW curves may coincide because the grain 
size is much smaller than the flux trapping depth^^ inside the sample. 


Figure 4(b) shows the resistive transition from the normal state to the superconducting 


state in the of the Nb 75 Zr 25 sample in various applied helds. The superconducting transi¬ 
tion width is less than O.IK in zero held which is quite unusual for the sample with such 
composition variation shown earlier. This further conhrms that though there is substantial 
amount of disorder in the sample, it behaves like a bulk superconducting sample of single 
composition. The width of the transition as a function of held is shown in the inset to hg. 


4(b) The observed width can be described by the equation. 


ATc = ATc{0) + FC 


( 1 ) 


with c = 1.82. The value of the exponent is not 2/3 as calculated theoretically for a giant hux 
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(a) 


(b) 


FIG. 5: (a) Resitivity as a function of inverse temperature at various applied fields. The solid 
lines are fit to the non-Arrhenius relation mentioned in equation 01 (b) Resistivity as a function 
of field at various temperatures. The solid lines are fit to the non-Arrhenius relation mentioned in 
equation [5l 


creep model or superconducting glass.- However, it should be noted that the quantitative 
dependence of the transition width on the applied held depends on the nature of activation 
energy.- The quantitative predictions of the theory may actually depend on the huxon 
jumping length.^ A deviation from the 2/3 behaviour has been observed in other systems 
as well.— >21 It is the central aim of this work to hnd out the distance over which the hux 
moves before a hnite resistance appears across the transition. It is therefore important to 
know the exact form of the activation energy as we see next. 


Figure 5(a) shows resistivity as a function of inverse temperature. The resistivity cannot 


be explained by a simple arrhenius relation where the activation energy is independent of 
temperature. Rather, it takes a non-arrhenius form which can be expressed as,— 




( 2 ) 


where 

U{H,T) = Uo{H){l-TlT,r (3) 

Substituting equation [3] in equation [2l we get the temperature dependent resistivity at a 
particular held as, 

p{T) = p,eM-U(T=o){H){l-T/TXIkBT] (4) 
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The best fit is obtained for a = 4 with a regression better than 0.996. Figure 5(b) shows 
the resistive transition as a function of magnetic held at various temperatures. The magnetic 
held dependent isothermal resistivity can be expressed with a similar non-arrhenius form as, 


p(H) = poexp[-t/(H=o)(T)(l - H/Hc2f/kBT] 


( 5 ) 


with the exponent /3 once again = 4. The activation energies U{H, T) determined from both 
temperature dependent resistivity in constant magnetic held and isothermal held dependent 
resistivity match quite closely as can be seen from table HI The prefactor Uq is quite large 
and is not the actual activation energy. Rather, it is the tangential slope of the arrhenius 
plot and is thus termed as the apparent activation energy.— 


Figure 6(a) shows the temperature dependent activation energy at various helds. The 
T(eV) line separates the hux creep and the thermally activated hux how regions. When 
U{T) << ksT, the thermally activated hux how is the dominant phenomenon responsible 
for hnite resistance. The variation of the prefactor Uq{H) in equation |3] is shown in hgure 


6 (b) The held dependence of Uo can be explained as a power law with a very large exponent. 


We obtain a value of -2.82 for the exponent which is considerably larger than what has been 
seen in other systems.— The signihcance of such a large value is not clear at present. 

The non-Arrhenius shape of the resistive transition is similar to one of the universal char¬ 
acteristics exhibited by glass-forming liquids, which is the non-Arrhenius viscosity across the 
glass transition.— >2^ The Arrhenius temperature dependence of viscosity arises when there 
is a temperature independent barrier which has to be overcome (known as the activation 
energy) for the flow of molecules to occur. However there are only a few exceptions like 


TABLE I: Activation energy determined from temperature dependent resistivity (R-T) in various 
fields and field dependent resistivity (R-H) at various temperatures. 


T(K) 

PoR(T) 

[/i?_r(meV) 

[/^-^(meV) 

10 

1 

1.465 

1.465 

9.5 

1.5 

3.374 

3.671 

9 

2 

8.172 

7.211 

8.5 

3 

0.0626 

0.0626 

8 

3.5 

0.468 

0.4775 

7.5 

4 

1.128 

1.117 
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(a) (b) 

FIG. 6: (a) Activation energy as a function of temperature at various fields as determined from 
the fit of equation 0] to the experimental data, (b) The prefactor Uq in equation 0] as a function of 
field. 

Si02 and Ge02 which exhibit an Arrhenius viscosity and in majority of cases, the activa¬ 
tion energy increases upon cooling.— A non-Arrhenius viscosity can also arise in case of a 
“landscape” of potential energy in the conhguration space.— 

The shoving modef^ describes the non-Arrhenius temperature dependent viscosity of the 
glass forming liquids. We hrst briefly describe this model and then show how this model, 
which was originally proposed to explain the glass transition in molecular liquids, can be 
actually applied to the case of vortex matter to explain our observed results. We reproduce 
here some of the arguments by Dyre et. al^ for the sake of continuity. 

The starting point of the shoving model is to view viscous liquids as ‘solids which flow’. 
Such a solid flows under force by sudden, rare, and localized molecular rearrangements and 
extra volume is needed for the flow event to occur. The work done in creating this extra 
volume is the activation energy and this is how the name of the model arose: in order to 
rearrange, the molecules must shove aside the surrounding molecules.— The elastic constant 
determining the shoving work is the shear modulus. The measure of how much the activation 
energy changes with temperature is the ‘index’, I = -dln[//dlnT and ranges between 2 and 
6 for most of the glass forming liquids. In our case, the value of the index is 4 if we take the 
logarithmic derivative of the activation energy with respect to the reduced temperature as 
can be seen from equation 01 

The temperature dependent activation energy in the shoving model is given by 
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U(T) = V^MT) 


( 6 ) 


where Vcorr is a correlation volume and G is the shear modulus. 

The basic assumptions behind Eq. [6] are as follows,— 

1. The activation energy is (mainly) elastic energy. 

2. This elastic energy is located in (mainly) the surroundings of the rearranging 
molecules. 

3. The elastic energy is (mainly) shear elastic energy, i.e., not associated with any density 
change. 

The model which was proposed for the molecular flow in viscous liquids (disordered solids) 
can be applied to the vortex matter if we replace the molecules with flux-lines. The flux- 
lines form a solid (lattice) under repulsion and thus have to shove aside the surrounding 
flux-lines if a flow has to take place. It is thus easy to see that the elastic energy is mostly 
located in the surroundings of the rearranging flux-lines. The elastic energy is also mainly 
the shear elastic energy as no density change is involved during the flow process because 
the density of flux lines is solely dependent on the applied held. Moreover, it costs a lot 
more energy to tilt the huxons as we shall see later when we calculate the elastic constants 
of the hux-line lattice. Therefore the only drastic assumption needed to apply the shoving 
model to the how of hux-lines is that the activation energy is mainly elastic energy. The 
Kramer model which explains the scaling laws of hux pinning implies this assumption while 
relating the strength of pinning with the appropriate elastic constant.— In this work we 
build upon those assumptions and relate them with the molecular how of disordered solids 
as described by the shoving model. We view the irreversible to reversible transition (or the 
resistive transition) in vortex matter as an event which can be described as the viscous how of 
disordered solid. In that respect, the various pictures of glass-like behaviour,->^ hux creep,-*^ 
or the melting of hux line lattice^*^ can be treated as equivalent, except for the shape of 
the irreversibility line, with the shear modulus of the vortex lattice as the key ingredient 
behind all these observed phenomenon. The apparent similarities (and diherences) between 
the various interpretations of the resistive transition in vortex matter have been discussed 
in details by Brandt.— 
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For the following part of this article, the focus will be on equation [6l The activation 
energy is already obtained from resistivity measurements. Our task is now to obtain the 
shear modulus of vortex matter for our sample and identify a volume which can be used as 
the correlation volume. We use the magnetization measurements for this purpose. 


C. Magnetization measurements and determination of elastic constants of vortex 
matter 


Figure [7] shows the isothermal magnetization as a function of held of the Nb 75 Zr 25 sample 
at T=7K. The insets show the various superconducting parameters obtained from this curve. 
The mixed state is paramagnetic at high helds and substantial positive magnetization can 
be observed even in the superconducting state. The upper critical held Hc 2 is determined 
as the held value at which the magnetization deviates from the linear M-H curve, passing 
through origin, of the paramagnetic phase in the normal state. The irreversibility held is 
estimated from the opening of hysteresis between the held-increasing and held-decreasing 
curves. 

The equilibrium magnetization is estimated as. 


where M 7 is the magnetization of the held-increasing cycle and M is the magnetization 
on the held-decreasing cycle. 

The critical current density for the rectangular sample (in SI units) is estimated as,— 


MH) 


AM{H) 
02(1 — 02/301) 


( 8 ) 


where AM is the diherence in magnetization between the held-decreasing and held- 
increasing curves at a particular held value and ai and a 2 are sample dimensions perpendic¬ 
ular to the direction of applied held, with oi > 02 . 

The thermodynamic critical held He is calculated from the equilibrium magnetization as 


rHc2 


Ht, = 2 


M,,{H)dH. (9) 

Jo 

The critical helds and the irreversibility line are plotted in hgure [HI It should be noted 
that the irreversibility line does not follow the 2/3 dependence which is seen in case of the 
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FIG. 7: Isothermal magnetization curve at 7K. The insets show the determination of various 
physical parameters from the magnetization measurement, (a) Determination of irreversibility 
field, (b) Determination of upper critical field Hc 2 - (c) The equilibrium magnetization as a 
function of field as determined from equation [7l (d) The critical current density as a function of 
field as determined from equation [8l 

vortex glasa^ or the giant flux creep.* Rather it is seen that (1 — T/Tc) oc which 

does not match with any of the models to explain the resistive transition in vortex matter. 
However, as we have pointed out earlier, the shape of the irreversibility curve is not the 
focus of our present work and we are mainly concerned about the distance over which the 
flux lines move when a hnite resistance appears across the transition. 

Apart from the critical helds, the other useful quantities are the fundamental length scales, 
the coherence length ^ and the penetration depth A, which will be useful to understand the 
pinning properties when compared with the microstructure discussed earlier. The coherence 
length is estimated from the upper critical held as, 

.2 ^ 4>0 

2tiijloHc2 

The penetration depth is estimated from its relation to the lower critical held as. 



Hci = 


47r/roA2 ^ ’ 

where k is the GL parameter given as X/^. The penetration depth, coherence length and 
the GL parameter are shown in hgure [H 
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FIG. 8: The field-temperature phase diagram of Nb 75 Zr 25 which shows the upper critical field and 
the irreversibility line. Inset shows the lower critical field and the thermodynamic critical field as 
a function of temperature. 


The interesting point to be noted here is that the coherence length is quite large as 
compared to the typical defect sizes shown in figure [3l The penetration depth also being 
quite large, allows the field to penetrate over a large distance around the vortex core and 
the held may actually span a lot of defects before decaying. It is therefore more likely 
to have more than one pinning mechanism as we see later during the description of the 
pinning properties. The large coherence length could also be the reason behind the sharp 
transition observed in bulk measurements. The coherence length for the composition which 
has the largest Tq could turn the surrounding ‘normal’ regions into superconducting regions 
through proximity effect, thereby making the composition inhomogeneity quite ineffective 


in broadening the transition. Figure 9(b) shows the variation of the GL parameter as a 
function of temperature. There is a sharp decrease in the value of k near 8K whose reason 
is not clear at present. 

The disorder profile in our sample also gives rise to the possibility of collective pinning 
of the flux-lines. If the flux-lines are rigid and placed in a perfectly periodic arrangement, 
the pinning would be quite inefficient due to the random nature of the underlying pinning 
centres, which are in the form of composition variation, dislocations and amorphous regions 
of the dendritic arms (see figures [2] and [3]) . It is therefore favourable for the individual flux- 
lines to lower their energy by passing through the nearest pinning sites and deviate from the 
ideal periodic arrangement. This distortion (deformation) of the flux-lines however increases 
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FIG. 9: (a) The characteristic length scales, penetration depth (A) and coherence length {^) as a 
function of temperature and (b) The GL parameter (k) as function of temperature for Nb 75 Zr 25 . 

the elastic energy and the eqnilibrium flux-line conhguration will be that arrangement which 
minimizes the sum of these two energies. This is the central theme of the theory of collective 
pinning by Larkin and Ovchinnikov.— The distortion in the flux-line lattice can be then 
understood in terms of a certain correlation radius and correlation length within which the 
lattice remains reasonably undistorted. The elastic energy consists mainly of the shear and 
tilt experienced by the flux line lattice due to this distortion. The shear modulus Cqg is 
given by,— 

where the reduced held b = B/Bc 2 - The tilt modulus 044 is given by, 

044 = BHa (13) 

Figure fill] shows the shear modulus and the tilt modulus as a function of reduced held 
at various temperatures. It can be seen that the tilt modulus is nearly 3 to 4 orders of 
magnitude larger than the shear modulus for most of the held range. This implies that it is 
much easier to shear the hux lattice than to tilt it. This observation is important because 
the shoving model assumes that the how occurs mostly by shearing the solid, which is the 
basis of equation [HI 

We shall hrst check the applicability of the collective pinning theory to the present 
Nb 75 Zr 25 sample by studying the pinning properties. Once the applicability is justihed, 
the correlation lengths obtained from the collective pinning theory along with the elastic 
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(a) (b) 

FIG. 10: Elastic constants of vortex matter as a function of reduced field at various temperatures, 
(a) Shear modulus and (b) tilt modulus. 

constants of the flux line lattice will be used to solve equation |6] in terms of the experimen¬ 
tally observed parameters. 


D. Pinning properties and the applicability of collective pinning theory 

To check if the pinning properties can be explained within the framework of single vortex 
pinning or collective pinning, the estimate of the depairing current is quite useful. The 
depairing current is given by,— 


Jn 




c 


4H, 


c 


\/6/ioA \/6A 


(14) 


The value of the depairing current at 7K (estimated from figures [H] and 9(a)) turns out 
to be 4.98 x 10^^ Am“^, with the highest value of Jc being of the order of 10® Am“^. 

The correlation length for single vortex pinning is,— 


T sv c 



1/2 


which at 7K and 4T is 3.39 x 10 '^m. 

The field region below which the single vortex pinning is applicable is,— 


(15) 


Bsv ~ ho 



(16) 


19 









FIG. 11: (a) Normalized pinning force density as a function of reduced field at various temperatures 
and (b) Kramer plot . 


which is 2.9 x 10“^T at 7K. This field value is much smaller than the value of lower 
critical field. Moreover, the ratio of the critical (or the depinning) current density and the 
depairing current density is quite small, which allows us to discuss the pinning properties in 
terms of the weak collective pinning theory.— 


Figure 11(a) shows the normalized pinning force density as a function of reduced field. 


Two distinct temperature regimes can be identified in fig. 11(a) where the shape of the 
pinning curves are quite different. At temperatures below nearly 8K, the curve is quite 
broad and appears to be a combination of two pinning laws. At temperatures above 8K, 
the peak in the curve shifts to much lower field values. The distinction of the temperature 
regimes also coincide with the sharp change in the GL parameter as a function of temperature 


(see fig. 9(b)). Similar shift in the peak of the pinning curves with temperature has been 
observed in case of polycrystalline MgCNia with graphite nanoprecipitates.— The Kramer 
plot is shown in figure 11(b) The kramer plot as a function of B) is piecewise 

linear with two major straight line portions which indicates that there are at least two 
dominant pinning mechanisms present in our sample. The normalized pinning force density 
can be scaled with the combination of two pinning laws as. 


f(p) = + C2ifyi - ( 17 ) 

Pmax 

Figure [12] shows the normalized pinning force density at two temperatures, one far below 
Tc and one at a higher temperature near Tc- The pinning curves can be fitted quite well 
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(a) (b) 

FIG. 12: The normalized pinning force as a function of reduced field at two different temperatures 
along with the fit to equation [T71 The total of the two scaling laws and the individual components 
are plotted along with the experimental data. 


with equation [T71 Such a combination has been used earlier to explain the pinning curve in 
Chevrel phase superconductors.— It was shown that two pinning mechanisms were possible 
within the framework of collective pinning depending upon the relationship between the grain 
size and the correlation radius.— The correlation radius in the collective pinning theory is 
given by,— 


Rc 




BJ, 


c 


(18) 


and the correlation length is 



(19) 


where, 


ttfu = 1.075 X 



( 20 ) 


is the flux line lattice parameter for a triangular Abrikosov lattice.— 

The variation of the correlation radius Rc and the correlation length Lc as a function 
of reduced field is shown in hgure [13] at three representative temperatures. The correlation 
radius is in a plane perpendicular to the direction of applied field and the correlation length 
is along the flux line. As we have seen earlier (fig. [TU]) . the tilt modulus is much larger 
than the shear modulus and hence the pinning properties could be quite insensitive to the 
correlation length and mainly depend on the relative changes of the correlation radius. The 
correlation radius remains fairly constant till a reduced held of about 0.7 at all temperatures. 
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(a) (b) 

FIG. 13: (a) The correlation radius and (b) the correlation length as a function of reduced field at 
three representative temperatures. 


At slightly higher helds, there is a drastic increase by almost an order of magnitude before 
Rc tends to drop down to zero at Bc 2 - The observed two pinning regimes can be explained 
with this rapid change in Rc as we discuss next. 

Figure [12] shows the normalized pinning force density at 4K and 9K each htted with 


equation [13 The normalized pinning force at 4K (£g. 12(a)) clearly shows the presence of 
two pinning mechanisms, one with the pi ~ 0.5 and qi ~ 2, which is commonly observed for 
grain boundary pinning or surface pinning^ with a peak at b ~ 0.2. The second component 
of equation Ogives quite high values of p 2 and q 2 - The htting is done till b ~ 0.8 by following 
the conventionally accepted procedure as the inhomogeneities in the sample signihcantly 
affect the tails of the pinning curves.— Such high values of the exponents have been reported 
in other systems as well and are thought to arise due to inhomogeneity in the sample.— 
Treating the exponents p 2 and q 2 as free running parameters in the htting procedure yields 
the best £t as seen in case of Chevrel phase compounds^ and Ti-V alloys^ for example. 


At higher temperatures (£g. 12(b)) the peak of the pinning curve distinctly shifts to¬ 
wards b ~ 0.2 and the surface (or the grain boundary) pinning appears to be the dominant 
mechanism of hux pinning. 

The answer to why the surface pinning becomes more dominant at higher temperatures 
can be obtained from the behaviour of Rc as a function of temperature and its relation with 


the microstructure of the sample. From hgure 13(a) it can be seen that Rc changes by almost 
an order of magnitude (especially at higher fields) as a function of temperature. At lower 
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TABLE II: The fitting parameters in equation [T7] along with the correlation radii at the peak 
reduced-field values of the pinning curves. 


r(K) 

Cl 

Low field peak 

Pi Q.1 bfYiax 

Ac (pm) 

C2 

High field peak 

P2 0.2 braax 

Ac (pm) 

4 

3.453 

0.55 

2.05 

0.21 

1.21 

230.5 

4.99 

2.27 

0.539 

0.85 

5 

3.281 

0.52 

1.95 

0.211 

1.29 

183.2 

4.95 

4.15 

0.543 

0.917 

6 

3.414 

0.52 

2.05 

0.206 

1.424 

231.4 

5 

4.55 

0.524 

1.064 

7 

3.268 

0.48 

2.05 

0.189 

1.59 

240.8 

4.99 

4.80 

0.509 

1.241 

7.5 

3.275 

0.48 

2 

0.193 

1.761 

284.5 

5.05 

5.50 

0.478 

1.462 

8 

3.397 

0.48 

2.05 

0.191 

1.932 

133.0 

4.99 

5.05 

0.496 

1.644 

8.5 

3.401 

0.48 

2.05 

0.188 

2.24 

141.6 

4.99 

5.06 

0.496 

1.93 

9 

2.949 

0.4 

2.05 

0.162 

2.856 

157.6 

5 

5.30 

0.483 

2.575 


temperatures (~ 5K and below), Rq is nearly 1 /rm. With this correlation radius, the defects 
and composition inhomogeneity within the dendritic arms (see £g. ED can offer effective 
pinning centres like the intragranular defects in case of Chevrel phase superconductors.— 
The dendritic arms which have a higher Zr concentration offer additional surfaces similar 
to grain boundaries in the average sample matrix. These additional surfaces could be the 
cause behind the usual pi = 0.5 and gi = 2 pinning law even when there are no well 
defined grains in our sample as we have seen during the discussion on microstructure earlier. 
The inhomogeneities within the dendritic arms could be the probable reason for the higher 
values of exponents p 2 and q 2 in equation [T71 At higher temperatures (and higher helds), the 


correlation radius grows by nearly 8 to 10 times (see fig. 13(a)) and is probably much larger 
than the internal structure of the dendritic arms. The inhomogeneities and defects within 
the dendritic arms are thus averaged out which makes the pinning over smaller length scales 
quite ineffective. This explains the shift of the peak in the pinning curve to lower helds and 
the dominance of the surface pinning at higher temperatures. 

The values of the htting parameters along with the correlation radii at the peak reduced- 
held values h^ax of the pinning curves are given in table [TTl It is thus justified that the theory 
of collective pinning can be applied to our sample and we now proceed to use this framework 
in equation [6] to answer the central question posed in the introduction. 
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E. Estimating the distance of flux movement during resistive transition 

We now compare the estimates of the displacement of the flux-line lattice by considering 
two processes of flux movement. One is the process of flux creep and the other is the process 
of the flow of viscous liquid as described in the shoving model. 

The pinning properties have shown that the influence of quenched disorder can be ex¬ 
plained in terms of collective pinning theory. To estimate the distance of flux movement 
when a finite resistance appears across the resistive transition, it is important to first know 
the volume of the flux bundle which is involved in the displacement. A collectively pinned 
object is classified as large-hundle when Rq > A, which seems to be the case in the present 
situation. 

The total elastic energy involved during the creep process can be estimated within the 
idea that a bundle of bundles, a superbundle with dimensions R\\ (parallel to the jump 
direction), R± = Re, (transverse to the held and the jump direction), and (along the 
held direction), will constitute the elementary unit.— (Details can be seen in figure 19 of 
Blatter et al.— ) By taking into account the compression of the neighbouring volume of the 
flux bundle during the jump process, one obtains a scaling relation,— 

R\\ = —Rl (21) 

afii 

When the tilt and compression energies are similar, R\\ ~ which are the longest 
dimensions of the superbundle.— 

For example, at 8 K and 2.9T in the case of the present sample, R± ~ Rc = 6.042 fim and 
R\\ ~ = 26.59 fim. Similar order of magnitudes for these characteristic length scales have 

been actually experimentally estimated on another conventional low-Tc material 2 if-NbSe 2 
by using neutron scattering.— 

The activation energy of a superbundle in the creep process near is given byi^ 

Ul ~ ( 22 ) 

where u is the displacement, which we term as Ax from here onwards. 

We wish to emphasize here that these are only order of magnitude estimates and the 
numerical factors have been dropped out from the calculations.— For estimating the dis¬ 
placement of the flux bundle during the creep process, the activation energy estimated from 
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resistivity measurements and the shear modulus estimated from the magnetization measure¬ 
ments is used. The typical length scales mentioned in these calculations and the associated 
displacement at various held and temperature values, when a hnite resistance appears just 
above the noise hoor across the resistive transition are given in table IIIIl 

To solve equation [H] we need the activation energy U, the correlation volume Vcorr and 
the shear modulus Cee- The Vcorr should not be confused with the correlation volume Vc as 
mentioned in the theory of collective pinning.— 

The correlation volume during shoving is given by the shoving model as,— 

~ {AVf/V (23) 

Where Vcorr = U/Cqq from the shoving model (equation [6]). (we have dropped the 2/3 
factor arising due to spherical symmetry to get only an order of magnitude estimate) V is 
volume indulged in shoving and AV is the volume change during shoving. We have taken 
U to be same as which is obtained earlier from £g. O The task is now to identify the 
volume V which is involved in shoving. The natural choice of such a volume within the 
collective pinning theory would be a parallelepiped of sides Rc (or R±) and length 
thereby giving a volume of R^L^ which would move when an external force is applied to it. 
In the shoving model, the shear energy is the only elastic energy related to the activation 
energy and thus the compression of the flux-line lattice is neglected. Therefore both the 
sides of the parallelepiped, R± and R\\ are taken to be same. 

To know about the change in volume during the shoving process, (i.e. the distance over 
which flux lattice moves) we have chosen those temperature and held values after which a 
hnite resistance appears just beyond the noise hoor during resistivity measurements. As the 
value of Jc enters into the calculation of Rc, the hrst non-zero value of Jc beyond noise 
hoor during magnetization measurements is used where, strictly speaking, the sample is not 
in the resistive state. Moreover, the appearance of hnite resistance is also a function of the 
sensing current through the sample. Therefore, the calculation is only an order of magnitude 
estimate of the actual phenomenon and should be treated likewise. 

At T = 9K and HqH = 1.9T, we have Rc ~ 7.541 pm and ^ 29.82 pm. This gives, 

V = R^L^ ^ 16.95 X 10“^® (24) 

At these temperature and held values, U = 2.47 x 10“^^J and Cqq = 30.55 TAm“^ (or 
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Jm ^). This gives from equation [6l Vcorr = 8.08 x 10 m^. 

From equation [23] above, 

{Avy ~ Korr X 1/ ^ 1.37 X 10-V (25) 

This gives AV ~ 3.7x10“^® m^. The change in area during the shoving process would 
therefore be AV/L^ which is nearly 124.13x10“^® m^. From simple geometric considerations, 
when the flux-line lattice shoves the neighbouring area aside by a distance of Ax, the change 
in area AV /for a square of sides Rc turns out to be, 

A{Axf + AAxRc (26) 

This gives Ax ~ 4.11xl0-^°m. 

The displacement estimated from the superbundle creep picture turns out to be 
8.27xl0-^°m which is nearly double but of the same order of magnitude as estimated by 
using the shoving model. 

Table IIIII shows the comparative study of the flux-line lattice displacement within the 
shoving model^ (i.e. if the vortex matter is treated similar to the viscous flow of disordered 
solid) and the superbundle creep picture which is more conventionally used to explain the 
resistive transition in vortex matter— at various fields and temperatures at which the irre¬ 
versibility transition occurs. We see that in all cases, the estimates from both the pictures 
have the same order of magnitude. It should be noted that these are only order of magnitude 
estimates and the actual numerical factors are required (not available within theory, to the 
best of our knowledge) to compare these displacements with the flux-line lattice constant. 

The implications of our study are quite interesting. We have attempted to find the 
analogies between the flow of disordered solids or viscous liquids with the resistive transition 
in vortex matter. We have not assumed any model like the superconducting glass or lattice 
melting for explaining the shape of the irreversibility line. Our arguments are based purely 
on thermodynamic considerations of elastic energies and could thus be a special case of a 
more general phenomenon of pinning and depinning properties elastic manifolds in random 
media.— 
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TABLE III: The comparison between the displacement of flux line lattice estimated from the 
shoving model and within the superbundle creep picture at various temperature and field values 
along the irreversibility line. 


T 

(K) 

IXqH 

(T) 

U 

(10-21J) 

Cm 

(T.Am-i) 

Ri. 

(dm) 

R\\ 

(dm) 

(dm) 

Lc 

(mm) 

(nm) 

^^{shoving) 

(nm) 

^^(bundle) 

(nm) 

9.5 

1.3 

4.35 

18.28 

6.06 

24.1 

24.1 

2.33 

42.7 

0.78 

1.57 

9 

1.9 

2.47 

30.55 

7.54 

29.8 

29.8 

3.27 

35.4 

0.41 

0.83 

8.5 

2.3 

2.84 

80.84 

4.68 

18.8 

18.8 

1.51 

32.2 

0.34 

0.68 

8 

2.9 

2.54 

143.64 

6.04 

26.6 

26.6 

1.84 

28.7 

0.20 

0.39 

7.5 

3.6 

1.28 

118.11 

8.18 

37.2 

37.2 

3.42 

25.7 

0.13 

0.25 

7 

4 

3.10 

138.77 

7.35 

34.0 

34.0 

3.15 

24.4 

0.20 

0.38 

6 

5 

2.89 

185 

4.2 

19.6 

19.6 

1.95 

21.8 

0.22 

0.41 


IV. CONCLUSION 

In conclusion, we have studied the resistive transition of vortex matter in highly strained 
sample of Nb 75 Zr 25 . The detailed study of microstructure of the sample enabled us to apply 
the theory of collective pinning to understand the pinning properties of flux line lattice in 
our sample and determine the elastic constants of vortex matter. The non-Arrhenius shape 
of the resistive transition, both as a function of temperature and magnetic held, showed 
that the models of viscous how of disordered solids can be indeed applied to vortex matter. 
The arguments were based on purely elastic energy considerations instead of assuming any 
particular model or shape of the irreversibility line in the held-temperature phase space. Our 
results show that the viscous how of disordered solids and the hux-creep phenomenon in hard 
type-II superconductors could be the manifestation of same underlying physical principles. 
These studies should provide sufficient interesting inputs for further experiments which can 
image the actual how of hux line lattice and those theories which treat the hux how in 
type-II superconductors as special case of the general phenomenon of plastic depinning of 
driven systems. 
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